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For  the  past  four  years  spectroscopic  measurements 
and  modeling  activities  directed  to  developing  a more 
precise  understanding  of  laser  transmission  have  been 
in  progress  [l-4].  The  measurements  objectives  of  this 
work  has  been  to  develop  a spectroscopic  data  base  from 
which  atmospheric  absorption  line  strength,  width  and 
position  parameters  may  be  extracted.  The  long  term 
objective  has  been  to  use  this  data  base  to  develop 
modeling  algorithms  to  predict  atmospheric  molecular 
absorption  at  laser  frequencies.  It  has  long  been 
known  [2,5]  that  molecular  absorption  may  change  by 
orders  of  magnitude  in  a spectral  interval  of  less  than 
a wavenumber,  and  that  a precise  line  parameter  data 
base  is  required  for  satisfactory  predictive  modeling. 

Measurements  have  been  restricted  to  high  resolution 
(<  . 05  cm  ^ ) molecular  line  absorption  in  the  .9  vim, 

3-5  v>m»  and  10-11  vim  spectral  regions.  Most  recently, 
these  activities  have  been  restricted  to  the 
3. 0-5.0  vim  region  where  the  DF  laser  operates. 

The  investigation  reported  here  completes  the  measure- 
ments activities  directed  toward  developing  field  verified 
predictive  modeling  of  DF  laser  atmospheric  molecular 
absorption.  In  this  report,  some  181  HDO  line  profile 
measurement  sets  are  presented  for  the  path  length  and 
partial  pressure  combinations  required  to  determine  the 
air  broadened  HDO  line  parameters.  Also,  synthetic 
atmospheric  molecular  absorption  spectra  have  been 
generated  for  the  several  molecular  absorbers  (N^O,  CH^, 
HDO,  CO2)  known  to  be  important  in  the  DF  laser  region. 
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The  DF  transmission  algorithms  have  been  presented  in 
a companion  report  [4].  Comparisons  between  field 
data  [6,7]  and  predictive  modeling  developed  here  have 
been  made. 

During  the  course  of  the  DF  transmission  modeling 
effort,  it  became  apparent  that  the  line  of  sight  mole- 
cular content  should  depend  on  the  location  and  environ- 
mental conditions,  in  general,  and  that  techniques 
must  be  developed  to  monitor  molecular  content  during 
laser  transmission  measurements.  Since  HDO  had  been 
identified  as  the  most  important  line  absorber,  and  since 
its  natural  abundance  and  its  variability  with  local 
conditions  was  not  understood,  a Gas  Filter  Correlation 
Spectrometer  (GFCS)  was  designed,  constructed  and  delivered 
to  NRL  to  be  used  to  monitor  actual  line  of  sight  HDO 
content  [p].  Since  the  variability  of  HDO  in  the  real 
atmosphere  was  shown  by  NRL  to  be  quite  dramatic  and 
important,  a design  study  to  expand  the  GFCS  capability 
was  initiated.  This  study  is  included  in  Section  5 of 
this  report. 

A status  summary  of  DF  transmission  modeling  is 
included  in  Section  6. 
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2 

ATMOSPHERIC  TRANSMISSION  CALCULATIONS 

It  is  required  that  accurate  transmission  algorithms 
be  developed  for  all  DF  laser  lines  from  which  significant 
power  is  extracted  in  the  high  energy  laser  (HEL)  devices. 

For  this  reason,  synthetic  atmospheric  transmission  spectra 
have  been  calculated  using  the  SAI  developed  code  SYNSPC  l 4) 
to  establish  the  importance  of  each  absorber,  and  to  identify 
discrepancies  with  laboratory  or  NRL  field  data.  In  this 
section,  a calculation  of  composite  transmission  for  the 
Mid-Latitude  Summer  (MLS)  model  atmosphere  [9]  and  individual 
components  are  described.  The  components  considered  are 
CC>2 » N^O,  HDO , and  CH^.  In  all  cases,  the  N_,  continuum  is 
included  since  it  is  a non  variable  component  at  sea  level. 

2.1  CALCULATIONS  FOR  THE  MLS  MODEL  ATMOSPHERE 

Synthetic  transmission  spectra  for  the  Midiattitude  Sum- 
mer Model  Atmosphere  were  calculated  for  the  nominal  5 km 
path  used  by  NRL  for  their  long  path  measurements.  Plots 
using  the  SAI  SYNSPC  code  are  given  in  Figures  la  - lq. 

The  SYNSPC  code  has  been  described  in  the  companion  report 
l 4 ] • Since  the  motivation  for  these  calculations  is  identi- 
fication of  absorption  features  and  determination  of  trans- 
mission magnitudes,  infinite  resolution  is  assumed.  The 
nominal  concentrations  for  mixed  gases  are  used,  as  indicated 
by  the  information  format  above  each  plot.  The  .03%  value  for 
the  HDO/H^O  ratio  is  used,  although  tentative  NRL  data  indi- 
cates that  this  ratio  is  variable,  and  sensitive  to  meteoro- 
logical conditions.  The  N.,  continuum  is  that  measured  by 
Burch  (10),  In  addition  to  CO^  line  absorption,  it  is  well  known 
that  COj  line  wing  absorption  is  sub-Lorentzian , and  that  the 
use  of  Lorentz  shapes  under  predicts  transmission  most  impor- 
tantly at  wavelengths  shorter  than  the  CC^  band  head.  In 
these  calculations,  Lorentz  shapes  have  been  used. 
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FIGURE  11.  FIGURE  1 CONTINUED 
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2.2  C02  TRANSMISSION  SPECTRA 

The  CO^  component  of  the  atmospheric  transmission  shown 
in  Figure  I has  been  calculated  spearately.  These  results 
are  shown  in  Figures  2a  - 2n.  N2  continuum  is  included. 

The  C02  continuum  absoprtion  caused  by  the  Lorentz 
wing  is  evident.  Also,  the  spectral  region  between  ~2. 390  cm 
and  —2406  cm  ^ is  seen  to  have  a small  computational  error. 
This  error  causes  the  unrealistic  appearance  of  the  line  ab- 
sorption superimposed  on  the  continuum  and  wing  absorption. 
The  cause  of  this  effect  has  been  determined,  and  corrections 
to  the  code  will  be  made. 
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FIGURE  2k.  FIGURE  2 CONTINUED 
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2.3  N20  TRANSMISSION  SPECTRA 

Transmission  plots  of  the  N20  component  of  the  atmospheric 
transmission  shown  in  Figure  1 are  shown  in  Figures  3a  - 3k. 

As  was  the  case  for  the  C02  calculations,  the  Burch  N2  contin- 
uum is  included  in  the  N20  calculations.  Several  panels  have 
been  omitted  since  there  is  no  detectable  absorption  by  N20 
in  the  region. 
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FIGURE  3g . FIGURE  3 CONTINUED 


FIGURE  3i . FIGURE  3 CONTINUED 
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2.4  WATER  VAPOR  TRANSMISSION  SPECTRA 

The  water  vapor  component  (plus  the  N2  continuum)  of 
the  atmospheric  transmission  plots  shown  in  Figure  1 are 
given  in  Figures  4a  - 4q.  Three  water  vapor  components  are 
included:  (1)  H^O  lines,  (2)  HDO  lines,  and  (3)  water  vapor 

continuum.  The  water  continuum  suggested  by  Burch  is  used, 
although  more  recent  data  soon  will  be  available. 

Inclusion  of  the  N^  continuum  here  provides  a clear 
description  of  how  the  absorption  components  are  modeled 
into  the  SYNSPC  code  (and  typically  in  all  such  codes) . 

Figure  4a  shows  the  sharp  cut-on  of  the  N~  continuum  at 

-1  *■ 

—2389  cm  , and  its  merging  with  the  HDO  lines  at  shorter 

wavelengths.  This  artificial  "step-function"  modeling  of  the 

N2  cut-on  is  simpler  computationally,  and  it  does  not  affect 

the  transmission  at  longer  wavelengths,  since  the  C02  in  the 

atmosphere  is  so  strongly  absorbing. 


FIGURE  4a.  ATMOSPHERIC  TRANSMISSION  FOR  THE  WATER  VAPOR 
COMPONENT  FOR  THE  INDICATED  CONDITIONS 
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FIGURE  4b 


FIGURE  4d . FIGURE  4 CONTINUED 
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2.5  CH.  TRANSMISSION  SPECTRA 

4 

The  transmission  spectrum  through  CH^  above  is  given 
in  Figures  5a  - 5f . The  N2  continuum  is  outside  the  spectral 
range  of  these  plots.  CH^  is  the  only  molecule  for  which 
there  is  serious  discrepancy  between  the  synthetic  spectra 
and  data  [3].  The  current  plots  were  made  using  the  January 
1977  AFGL  data  compilation  which  includes  a set  of  CH^ 
parameters  which  is  an  improvement  over  those  available 
earlier.  Since  CH^  is  of  minor  importance  at  sea  level, 
laboratory  measurements  were  not  performed.  However,  the 
effect  of  CH4  on  laser  transmission  at  higher  altitudes 
or  for  long,  slant  paths  is  expected  to  be  significant. 

The  abrupt  termination  of  CH4  absorption  in  Figure  5f  is 
caused  by  a lack  of  entries  in  the  January  1977  ver  on 
of  the  AFGL  data  compilation. 
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HIGH  RESOLUTION  HDO  LINE  PROFILE  MEASUREMENTS 


A large  quantity  of  air  broadened  HDO  data  were  obtained 
at  several  total  pressures,  water  vapor  partial  pressures,  and 
path  lengths.  All  measurements  were  performed  on  the  3 meter 
focal  length  Ebert  spectrometer  located  at  the  Physics  Depart- 
ment of  The  University  of  Michigan.  Details  of  the  measure- 
ments apparatus  and  experimental  procedures  were  explained 
in  detail  in  an  earlier  report  [2],  and  therefore  they  will 
not  be  repeated  here.  In  the  following  paragraphs,  the 
measurements  plan  will  be  described  and  the  HDO  data  will 
be  discussed. 


3.1  HDO  MEASUREMENTS  PLAN 

A total  of  181  sets  of  data  were  obtained  on  HDO.  De- 
scription of  the  181  sets  are  included  in  the  Appendix. 

The  data  were  taken  to  fulfill  four  requirements. 

These  are  as  follows: 

(1)  Strength  and  width  determination  (S,y).  These 
are  the  basic  data  obtained  at  several  pressure 
and  path  length  conditions  from  which  the  strength 
and  width  values  are  obtained. 

(2)  Vacuum  absorption  cell  measurements  (EMPTY) . This 
information  is  required  for  determining  the  cell 
mirror  reflectivities. 

(3)  Low  pressure  data  (SLIT) . These  data  are  required 
for  determining  the  spectrometer  slit  function,  for 
use  in  FITLINES. 

(4)  HDO  calibration  data  (HDO  cal) . These  data  were 
obtained  to  verify  that  the  H_0  + HDO  + D2O  mix- 
ture had  equilibriated  with  themselves  ana  with 
the  cell  walls.  When  the  cell  was  filled  with 
different  amounts  of  the  D isotope,  scans  were 
made  to  verify  that  the  HDO  did  not  change  with 
time,  and  that  new  fills  were  existent  with  old. 
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(5)  Line  position  data  ( v ) . This  data  was  taken  with 
the  pulsed  DF  laser  operating,  to  obtain  HDO  line 
center  positions,  relative  to  the  DF  laser  lines. 

Data  was  taken  for  three  total  air  pressures:  4 PSI , 

8 PSI,  and  14.7  PSI.  A complete  summary  of  the  data 
is  given  in  Figures  6,  7,  and  8,  respectively,  for 
these  cases. 

A "run"  is  defined  as  a set  of  scans  for  which  the  cell 
fill  is  not  changed.  A "run"  is  an  individual  spectroscopic 
scan,  usually  limited  by  the  data  which  could  be  obtained  as 
a single  tape  cassette.  A complete  description  of  the  indi- 
vidual runs  and  scans  is  given  in  the  Appendix. 

3.2  HDO  MEASUREMENTS 

Line  profile  measurements  have  been  performed  throughout 
the  entire  spectral  region  that  it  is  a significant  absorber 
of  DF  laser  radiation.  The  entire  data  accumulation,  for 
strength  and  width  as  well  as  for  diagnostic  requirements 
are  summarized  in  the  Appendix. 

Reproduction  of  all  the  data  in  the  form  of  transmission 
plots  would  be  costly  and  extremely  voluminous.  Therefore, 
a smaller  set  of  data  has  been  selected  for  plotting  in  the 
SAI  format.  Only  the  spectra  of  samples  having  a total  pressure 
of  ~1  atm  have  been  plotted.  This  choice  is  made  since  this 
condition  is  representative  of  the  sea  level  propagation  alti- 
tude most  important  to  the  Navy.  Theoretical  plots  have  been 
calculated  for  the  same  path  length  and  pressure  conditions 
as  the  data  measured.  Strength  and  width  parameters  for  the 
theoretical  plots  were  obtained  from  the  AFGL  data  tape, 

January  1977  version.  These  plots  are  reproduced  on  the 
same  page  as  measured  data  for  easy  comparison.  The  agree- 
ment is  quite  good  relative  to  earlier  comparisons.  This  is 
essentially  because  this  data  has  been  used  to  update  HDO 
parameter  modeling  in  the  AFGL  compilation  prior  to  its 
publication  [l4]. 
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The  data  selected  for  publication  and  the  calculated 
plots  are  shown  in  Figures  9-14.  Figures  9 and  10  repre- 
sent an  initial  cell  fill  of  9.41  torr  total  water  vapor,  for 
path  lengths  of  380  m and  20  m respectively.  Figures  11  and 
12  describe  initial  cell  fills  of  6.61  torr  total  water  vapor, 
for  20  m and  40  m path  length  respectively.  Figures  13  and 
14  are  for  initial  cell  fills  of  4.37  torr  total  water  vapor, 
for  380  m and  300  m path  length,  respectively.  Identifica- 
tion of  D20  lines  present  in  the  spectra  has  not  been 
attempted . 
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FIGURE  9d . FIGURE  9 CONTINUED 
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FIGURE  10a.  HDO  MEASUREMENTS  AND  COMPARISON  PLOTS  FOR 
9.41  TORR  TOTAL  WATER  VAPOR  FILL  PRESSURE 
20  M.  PATH 


FIGURE  10b.  FIGURE  10  CONTINUED 


FIGURE  10c.  FIGURE  10  CONTINUED 


FIGURE  lOe,  FIGURE  10  CONTINUED 


FIGURE  lOf.  FIGURE  10  CONTINUED 
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FIGURE  lOg,  FIGURE  10  CONTINUED 


FIGURE  lOi.  FIGURE  10  CONTINUED 


FIGURE  10 j.  FIGURE  10  CONTINUED 
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FIGURE  lOp.  FIGURE  10  CONTINUED 


FIGURE  11a.  HDO  MEASUREMENTS  AND  COMPARISON  PLOTS  FOR 
6.61  TORR  TOTAL  WATER  VAPOR  FILL  PRESSURE 
20  M.  PATH 


FIGURE  11c.  FIGURE  11  CONTINUED 
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FIGURE  12a.  HDO  MEASUREMENTS  AND  COMPARISON  PLOTS  FOR 
6.61  TORR  TOTAL  WATER  VAPOR  FILL  PRESSURE 
40  M.  PATH 


FIGURE  12b.  FIGURE  12  CONTINUED 


FIGURE  13a.  HDO  MEASUREMENTS  AND  COMPARISON  PLOTS  FOR 

4.37  TORR  WATER  VAPOR  FILL  PRESSURE,  380  M.  PATH 
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FIGURE  14a.  HDO  MEASUREMENTS  AND  COMPARISON  PLOTS  FOR  4 
TORR  WATER  VAPOR  FILL  PRESSURE,  300  M.  PATH 


TEMP  QO  i 
PATH  (KM) 


FIGURE  14b.  FIGURE  14  CONTINUED 


FIGURE  14d . FIGURE  14  CONTINUED 
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PLANNING  AND  ANALYSIS  SUPPORT 

Three  kinds  of  planning  and  analysis  support  have  been 
provided  to  NRL:  (1)  degraded  resolution  transmission  cal- 
culations in  support  of  NRL  field  Scanning  Michelson  Inter- 
ferometer (SMI)  measurements,  (2)  transmission  measurements 
to  support  future  NRL  laser  propagation  measurements, 

(3)  fixed  frequency  DF  laser  extinction  calculations.  These 
activities  are  discussed  in  the  remainder  of  this  section. 

4.1  SMI  CALCULATIONS  AND  COMPARISONS 

Calculations  for  comparison  with  NRL  SMI  measurements 
have  been  performed  during  the  NRL  field  measurements  pro- 
gram. SAI  has  used  the  SYNSPC  code  for  values  of  the 
variable  parameters  which  correspond  to  the  field  measure- 
ments conditions.  This  primarily  means  selecting  the 
spectral  resolution  (slit  width  in  the  plot  formats)  to 
match  the  SMI  resolution,  and  to  match  the  temperature 
total  pressure,  water  vapor  pressure,  and  path  length. 

HDO  quantity  and  concentrations  of  the  non  water  vapor 
mixed  gases  were  taken  to  be  the  values  used  in  the  AFGL 
data  compilation. 

A sample  plot  of  NRL  data  and  a comparison  plot  is 
shown  in  Figure  15  for  the  spectral  region  2480  cm  ^ - 
2800  cm  . NRL  data  is  on  the  lower  panel,  and  the  SAI 
calculation  for  the  indicated  conditions  is  shown  in  the 
upper  panel.  The  comparisons  are  quite  good,  as  viewed 
qualitatively.  Laboratory  compairsons  have  proven  that 
the  molecular  line  parameter  data  are  now  excellent,  so 
discrepancies  may  be  attributed  to  continuum  extinction, 
both  aerosol  and  molecular,  and  to  number  concentrations 
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FIGURE  15b,  FIGURE  15  CONTINUED 


in  the  atmosphere.  Recent  NRL  data  on  path  integrated  HDO 
measurements  have  shown  that  significantly  less  HDO  is  in 
the  atmosphere,  typically,  than  the  assumed  .03%  HDO/HHO 
ratio  115].  This  is  apparent  in  Figure  15b  only  with 
careful  scrutiny. 

4.2  TRANSMISSION  AT  CO  LASER  FREQUENCIES 

Transmission  in  the  spectral  region  encompassing  the 
lowest  lying  CO  laser  frequencies  has  been  calculated  for 
the  Midlatitude  Summer  Model  atmosphere,  for  three  CO 
molecular  concentrations.  Values  of  .225  parts  per  million 
(PPM)  .75  PPM  and  .075  PPM  have  been  used.  A synthetic 
spectrum  is  shown  in  Figure  16  for  the  lower  CO  bound. 
Tables  1,  2 and  3 gives  the  extinction  value  for  the  lower 
lying  CO  laser  lines,  and  the  absorption  coefficient  of 
each  molecular  absorber  that  is  significant  at  the  CO 
frequency. 

4.3  DF  EXTINCTION  CALCULATIONS 

DF  extinction  algorithms  have  been  determined  in 
detail  from  the  measured  HDO  data,  from  earlier  data  and 
using  values  in  the  AFGL  data  compilation.  These  have 
been  reported  earlier  14].  However,  during  the  course 
of  this  program,  preliminary  absorption  coefficient  values 
were  determined,  and  power  weighted  absorption 
coefficiencies  were  calculated  for  a measured  Baseline 
Demonstration  Laser  (BDL)  power  spectrum.  This  has 
been  compared  with  estimates  based  on  NRL  measured  extinc- 
tion coefficients.  The  results  are  shown  in  Table  4. 


FIGURE  16.  MONOCHROMATIC  ATMOSPHERIC  ABSORPTION  IN  THE 
CO  LASER  REGION 


Table  1. 

CO  Laser  Absorption  Coefficients 
Midlatitude  Summer  Sea  Level  with  0.075  ppm  CO 
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Table  2. 

CO  Laser  Absorption  Coefficients 
Midlatitude  Summer  Sea  Level  with  0.225  ppm  CO 
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2099.865 

255.9 

14.3 

1.5 

699.6 

2099.085 

69.  1 

2.  b 

2.0 

796.8 

PI  (10) 

2103.273 

215.2 

2.0 

1.6 

880.2 

PI  (9) 

2107.  926 

1369.0 

15.9 

2.  1 

939.6 

PI  (B) 

2111.596 

164.  1 

3.2 

0.2 

972.3 

PIP) 

2115.  632 

490.9 

B . 2 

1.2 

970.5 

PI  (5) 

2119.669 

35.2 

10.2 

125 

2.3 

92m  . 0 

TOTAL 

341.5 

5977.2 
1 6b7 . 8 

436.  4 
13*4 . 1 
5482.  4 

306.5 

1391.7 
326-  3 
870.  0 
326.9 

15860. 9 
10505.3 

300.6 

308.7 

1593. 8 
136.  B 
161.  1 

1288.8 

891.5 

1 0397.6 

102.7 
95. B 

58.2 
62.5 

366.  6 

37765.2 

1888.0 
177.1 
200.0 

1869.9 

5131.0 

387.9 

680.6 
1 183.0 
3019.  3 

700.  6 

2996.7 
971.3 

870.6 

1099.0 
2326.  6 
1 139.7 

1970.8 

971.7 


Table  3. 

CO  Laser  Absorption  Coefficients 
Midlatitude  Summer  Sea  Level  with  0.75  ppm  CO 


LASER 

SAI 

FREQUENCY 

ABSORPTION 

COEFFICI 
( 6-18- 

ENTS 

76) 

X 1000 

(KR-1) 

LINES 

(cr-1) 

820 

C02 

03 

CO 

P2(23) 

2011.419 

341.  4 

0.  1 

O 

• 

o 

0.0 

P2  (24) 

2016.025 

5977.;.0 

0.  1 

0.  0 

0.  1 

P2  (23) 

2020.  600 

1687.0 

0.  7 

0.  0 

0.5 

P2  (22) 

2025.  143 

435.4 

0.5 

0.0 

1.7 

P2  (21) 

2029.  655 

132.  6 

1.1 

0.  0 

1.3 

P2  (20) 

2034. 134 

5479.0 

3.3 

0.  0 

0.3 

P2  (19) 

2038.581 

291.  4 

15.0 

0.0 

0.2 

P2  (18) 

2042.996 

1173.0 

168.6 

0.  0 

0.3 

P2(17) 

2047.379 

259.2 

67.0 

0.0 

0.4 

P 2 (16) 

2051.729 

847.2 

22.5 

0.  0 

0.9 

P2  (15) 

2056.620 

69.  5 

25b.  9 

0.  1 

1.3 

P2 (14) 

2060. 332 

15340.  0 

18. b 

0.  1 

6.6 

P2 (13) 

2064. 583 

10380. 0 

112.  3 

0.2 

42.4 

P2(1  2) 

2068.802 

150.  7 

31.5 

0.  3 

393.5 

P2(  11) 

2072.987 

277.  9 

19. C 

0.7 

36.6 

P2  (10) 

2077.  139 

155.5 

1382.0 

0.7 

18.7 

P2  (9) 

2081.258 

119.  6 

9.8 

0.  9 

21.8 

P2  (8) 

2085.343 

137.  6 

19.8 

0.  9 

9.2 

P2  (7) 

2089.393 

1279.0 

6.  1 

1.4 

7.7 

P2  (6) 

2093.  410 

628.3 

260.0 

1.3 

6.5 

P 2 (5) 

2097. 393 

10390.0 

3.3 

2.  5 

6.3 

?2  (4) 

2101.392 

91.3 

7.9 

1.5 

6.5 

P2(3) 

2105.256 

88.8 

3.  1 

1.6 

7.1 

P2  (2) 

2108.768 

44.  4 

8.7 

1.8 

10.9 

P2  (1) 

2112.  549 

50.3 

6.  9 

0.  1 

17.2 

PI  (25) 

2 037.027 

327.  1 

22.0 

0.0 

58.2 

PI  (24) 

2041.668 

37720.0 

14.2 

0.  0 

103.  1 

PI  (23) 

2 046.278 

1829. 0 

20.0 

O.C 

129.8 

PI  (22) 

2050.856 

88.5 

32.7 

0.0 

185.9 

PI  (21) 

2055.402 

57.5 

63.2 

0.  1 

264.0 

PI  (20) 

2059.917 

1567.  0 

188.1 

0.  1 

365.7 

PI (19) 

2064.399 

4468.  0 

513.7 

0.2 

497.0 

?1  (18) 

2068.849 

149.  4 

39.5 

0.3 

662.5 

PI  (17) 

2073.267 

395.2 

25.5 

0.5 

865.0 

PI  (16) 

2077.652 

1 14.  6 

735.6 

1.0 

1106.0 

PI  (15) 

2082.005 

2588.0 

12.5 

0.7 

1377.0 

PI  (14) 

2086.325 

183.9 

10.9 

1.2 

1682.0 

PI  (13) 

-2090.611 

1842.  0 

2.3 

2.  1 

2001.0 

PI  (12) 

2094.865 

255.  9 

14.3 

1.5 

2332.0 

PI  (11) 

2099.085 

69.  1 

2.6 

2.  0 

2656.0 

Pi  (10) 

210  3.  27  3 

215.2 

2.0 

1.6 

2934.0 

PI  (9) 

2107.426 

136  9.0 

15.9 

2.  1 

3132.0 

PI  (8) 

2111.546 

164.  1 

3.2 

0.2 

3241.0 

PI  (7) 

2115.632 

490.9 

8.2 

1.2 

3235.0 

?1  (6) 

2119.  b84 

35.2 

10.2 

2. 3 

3080.0 
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TOTAL 

341.  5 

5977.2 

1688.2 

437.7 
135.  0 

5462- b 

306.7 

1341.9 
326-  6 
870. 
327. 

1 5865. 
10534. 
576.1 

334. 3 

1556.9 
152.  1 
1 o7 . 5 

1294-  3 
39b  . 1 
10402-0 

107.3 

100.  8 

65-8 

74.6 

407.3 

37837. 3 

1978.8 
307.  2 

384. 8 

2120.9 
5478.  9 

851.7 
1286.  1 

1957.2 

3978.2 

1378.0 

3847.4 

2603.7 

2729.8 
3152.  8 

4519. 0 

3408. 4 

3735.3 

3127.7 


vo  t/i  o:  O' 


TABLE  4 

BDL  WEIGHTED  a 

COMPARISON  OF  CURRENT  HI-TRAN  WITH  FIELD  MEASUREMENT 

P(H20)  -v  14  TORR 


Line 

% Power 

Calc  a 

Measured  a 

PI -5 

4.1 

0.070 

0.065 

PI -6 

6.1 

0.087 

0.074 

PI -7 

4.1 

0.064 

0.050 

PI -8 

14.8  - 

0.140 

0.152 

P2-5 

1.0 

0.030 

0.034 

PI -9 

0.5 

0.040 

0.038 

P2-6 

12.8  • 

0.059 

0.061 

P2-7 

11.3 

0.100 

0.114 

P2-8 

12.3 

0.035 

0.023 

P3-5 

0.8 

0.023 

0.019 

P2-9 

5.6 

0.041 

0.050 

P3-6 

4.3 

0.030 

0.032 

P2-10 

2.6 

0.068 

0.076 

P3-7 

12.8 

0.065 

0.068 

P3-8 

2.6 

0.051 

0.056 

P3-9 

3.6 

0.036 

0.033 

P3-10 

0.8 

0.051 

0.043 

Measured  a = 0.074 
Calc,  a = 0. 071 
Difference  = 4.3  % 


l 

l 
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5 

MODIFICATION  OF  THE  HDO  GFCS  TO  MEASURE 

ch4,  n2o,  c13o2,  and  c12o2 


NRL  measurements  of  path  integrated  HDO  vapor  have 
demonstrated  that  the  abundance  of  HDO  relative  to  HHO  is 
less  than  the  generally  accepted  value  of  .03%,  and  that 
it  varies  significantly  with  local  environmental  conditions  [15]. 
The  dramatic  effect  of  variations  in  molecular  content  on  DF 
line  absorption  naturally  leads  one  to  question  the  relia- 
bility of  accepted  concentrations  for  the  important  "uni- 
formly mixed"  gases  at  specific  sea  level  sites.  For  this 
reason,  a study  has  been  performed  to  determine  the  feasi- 
bility of  modifying  the  NRL  GFCS  to  measure  thevOther  gases 

important  to  DF  propagation.  These  gases  are:  CH.,  C330_, 

12  * * 

N20  and  C Oj • 

The  existing  GFCS  was  designed  and  built  under  Contract 
N00173-75-C-9627  and  delivered  to  NRL  in  October  1976.  The 
instrument  is  completely  described  elsewhere  [8]  . It  is 
an  infrared  instrument  designed  to  measure  the  atmospheric 
water  vapor  content  over  a multi-kilometer  path.  The  unit 
has  been  designed  to  operate  in  the  NRL  IMORL  field 
measurement  vans.  A source  unit  is  placed  in  the  trans- 
mitter van.  The  GFCS  unit  (Fig.  17) , a calibration  source 
unit,  and  an  electronics  unit  are  placed  in  the  receiver 
van  which  is  5 kilometers  away  from  the  transmitter  van. 

The  source  unit  provides  an  infrared  signal  which  is  sent 
through  the  5 kilometer  atmospheric  path  by  the  transmitter 
van  optics  to  the  receiver  van  optics.  The  GFCS  unit  then 
correlates  the  spectral  structure  of  the  received  signal 
with  the  spectral  structure  of  water  vapor  to  determine 
the  integrated  atmospheric  water  vapor  content. 
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The  suggested  modification  will  permit  the  measurement 

13 

of  N20,  CH4,  and  CC>2  over  the  same  5 km  path.  The  C C>2 
measurement  will  be  made  over  a short  path  of  approximately 
100  meters. 

5.1  ANALYTICAL  DESIGN  STUDY 

5.1.1  CONCEPT  OF  OPERATION 

The  GFCS  concept  was  explained  in  the  previous  report 
describing  the  instrument  and  should  be  referred  to  for 
more  detailed  explanation.  The  term  gas  filter  refers  to 
the  use  of  the  complex  infrared  absorption  spectra  of  a 
gas  as  an  infrared  filter.  Correlation  refers  to  the  spec- 
tral correlation  between  the  gas  filter  and  the  infrared 
absorption  spectra  of  the  same  gas  in  the  atmosphere.  The 
principle  of  operation  is  given  by  Equation  (1) . 


M 


T - T T 
ac  a c 


+ T T 
a c 


(1) 


where 

M is  the  instrument  signal  caused  by  the  gas  to 
be  measured 

T is  the  average  transmission  of  the  measured  gas 
a over  the  infrared  bandpass  of  the  GFCS 

T is  the  average  transmission  of  the  GFCS  sample 
cell  over  the  infrared  bandpass  of  the  GFCS 

Tgc  is  the  average  transmission  of  a combination  of 
u the  atmospheric  gas  transmission  and  the 

sample  cell  gas  transmission  over  the  infra- 
red bandpass  of  the  GFCS. 

Equation  (1)  expresses  the  measurement  signal  in  terms 
of  the  degree  to  which  it  modulates  the  source  signal.  The 
measurement  signal  M arises'  from  the  fact  that  the  average 
transmission  of  the  combination  of  atmosphere  and  sample 


r! 


i 


cell  (T  ) is  greater  than  the  product  of  individual  average 

transmissions  T and  T 
a c 


ac 


i Vc 


(2) 


where 


T 


ac 


Av 


T (v)T  (v)dv 

a C 


(3) 


Vc  * ^ [ TA,V’d'>  ” £ [ TCMdv  (4> 

JAv  JAv 

In  the  actual  GFCS,  the  source  signal  is  initially 
chopped  at  750  Hz  to  allow  the  signal  processing  electronics 
to  discriminate  between  it  and  other  signals.  This  beam  is 
alternately  passed  through  a sample  cell  to  obtain  the  com- 
bined cell  atmospheric  path  transmission  T , and  passed 
through  a reference  attenuation  path  with  a constant  trans- 
mission adjusted  to  the  average  sample  cell  transmission 
T_  = T to  obtain  the  product  T T . After  the  signal  is 
modulated  by  this  procedure  it  is  imaged  onto  a detector. 

The  modulation  M is  of  central  importance  because  the 
primary  limits  to  the  accuracy  of  the  measurements  are 
manifested  in  terms  of  the. modulation  of  the  source  signal 
rather  than  as  a constant  noise.  The  constant  detector 
noise  is  small,  less  than  the  amount  that  would  cause  a 
1*  measurement  error.  However,  errors  which  express  them- 
selves in  terms  of  the  modulation  limit  the  measurement 
accuracy.  The  purpose  of  the  analytical  design  studies  are 
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to  calculate  the  dependence  of  the  modulation  M on  the  various 
atmospheric  and  sample  cell  parameters:  spectral  region,  inter- 
fering species,  cell  concentration,  and  path  length. 

5.1.2  SPECTRAL  BAND  SELECTION 

A check  was  made  of  the  spectral  absorption  bands  of 
CH^ , H20  and  C02  which  might  be  suitable  for  a GFCS  instru- 
ment. The  major  bands  considered  are  listed  in  Table  5.  An 
asterisk  has  been  placed  next  to  the  selected  band. 

Accuracy  estimates  for  the  selected  bands  were  made 
based  on  line-by-line  calculations.  This  information  is 
given  in  Table  5. 

The  atmospheric  transmission  over  long  paths  for  the 
selected  bands  are  shown  in  Figures  5,  3,  and  18.  It  can 
be  seen  in  Figure  5 that  the  CH4  band  is  quite  strong,  but 
there  are  many  water  vapor  lines  in  the  same  region.  How- 
ever, our  calculations  predict  that  the  interference  will 
be  less  than  2%.  It  should  be  noted  in  Figures  3 and  18 
that  there  is  very  little  interference  in  the  N00  band  and 
virtually  none  in  the  CC>2  band. 

The  predicted  modulation  functions  for  these  same  bands 
as  a function  of  the  expected  concentration  is  given  in 
Figures  19  through  21.  These  can  be  compared  with  HDO  modu- 
lation function  given  in  Figure  22. 

In  addition  to  the  gases  included  in  the  detailed  analy- 
sis, it  is  proposed  to  measure  C^O,.  The  spectral  region 

1 -1 

selected  for  this  measurement  is  2240  - 2300  cm  . An  atmos- 
pheric transmission  path  in  this  spectral  region  is  shown  in 
Figure  23.  The  absorption  in  the  region  is  too  strong  to 
make  measurements  over  a sea  level  5 km  horizontal  path. 
Therefore,  it  is  being  proposed  that  a separate  portable 
source  be  provided  which  can  be  set  on  a tripod  at  a con- 
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TABLE  5 

GFC  BANDS  CONSIDERED  FOR' VARIOUS  MOLECULES 


Molecule 

pm 

Band  , 

cm 

Band  Strength 

i «.  “1  -2) 

(atm  cm 

Comments 

CH4 

1.665 

6005 

1.6 

Too  weak 

*3.312 

3019 

270.0  (a- 0.30) 

H-0  interference 

2860  - 2960  cm"1 

Not  bad. 

3.846 

2600 

3.07  (?) 

Probably  too  weak. 

' 

7.657 

1306 

204.0 

No  transmission  because 
of  H20  and  C02 

n2o 

17.007 

588.8 

35 

No  transmission 

HjO  and  C02 

8.561 

1168.1 

10 

Weak  interference 

: 

7.782 

1284.9 

230 

No  transmission 

4.497 

2223.8 

1600 

No  good  - badly  over- 
lapped with  C02 

*3.755 

2663.3 

40  (a  = 0.09) 

Looks  good,  may  be 
slightly  weak 

E 

2.873 

3480.8 

40 

No  transmission  HjO 

co2 

2.060 

4853.6 

0.217 

Looks  good  if  not  too 
strong 

2.080 

4807.7 

0.018 

Combination  band, 

4780-4804  fairly  clear 

• 

*1.606 

6227.924 

0.0137  (a  « 0.08) 

Clear,  no  interference 

1.576 

6347.  • 

0.0132 

. 

Clear,  no  interference 

• 

- 1 ’•  ■ 


133 


n 


TABLE  6 

ACCURACY  ESTIMATES  BASED  ON  LINE-BY-LINE 
CALCULATIONS  FOR  GFCS  GAS  CONCENTRATION 
MEASUREMENTS  OVER  A 10  KILOMETER  PATH 


Measured 

Gas  Measurement  Range  (ppm)  Accuracy 

C02  50  - 500  5% 

No  Interference 


N_0  0.05  - 5 10% 

No  Interference 


CH.  0.4-4  5% 

H20  Interference  < 2%  at  1.6  ppm 


i 
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venient  distance  from  the  receiver,  about  100  meters  for 
measurement.  Calculations  on  single  lines  in  this  region 
indicate  that  this  would  be  the  approximate  distance  re- 
quired. 


5.1.3  DETERMINATION  OF  CELL  CONCENTRATIONS 

The  modulation  function  M is  dependent  on  the  total 
amount  of  gas  in  the  cell,  the  temperature  and  the  pressure. 
Since  we  are  making  measurements  at  atmospheric  pressure, 
the  best  correlation  will  occur  with  lines  of  approximately 
the  same  width.  Hence,  the  first  choice  is  to  make  the 
cell  pressure  equal  to  atmospheric  pressure  and  maintain 
the  temperature  near  ambient. 

The  modulation  functions  plotted  in  Figures  21  - 23 
were  done  with  the  following  selected  cell  concentration: 

Molecule  Concentration 


CH4 

5 

atmos  - cm 

n2o 

3 

atmos  - cm 

C02 

1200  atmos  - cm 

The  estimated  concentration  required  for  C*"*02  is 
0.25  atmos  - cm. 

Before  the  final  cell  depth  is  selected,  the  modulation 
function  for  a family  of  cell  optical  depths  will  be 
calculated  similar  to  what  was  done  for  the  HDO. 

5.1.4  SELECTION  OF  MODULATION  SYSTEM 

As  was  pointed  out  in  the  introduction,  the  limits 
on  the  accuracy  of  measurement  are  set  by  the  modulation 
M rather  than  detector  noise.  In  the  case  of  the  HDO 
measurement,  the  nominal  value  of  M is  about  10%  whereas 
in  the  case  of  N20,  the  nominal  value  is  about  1%.  This 
means  that  in  order  to  achieve  the  same  ultimate  measurement 


ABSORPTION  THROUGH  A 10  KM  HORIZONTOL  PATH 
AT  12  KM  ALTITUDE 


accuracy,  the  minimum  value  of  M must  be  an  order  of  magnitude 
lower.  We  have  therefore  given  consideration  to  other  modula- 
tion schemes  in  addition  to  the  chopper  modulation  used  in 
the  current  GFCS. 

Two  alternate  modulation  schemes  were  considered.  The 
first  splits  the  incoming  radiation  with  a beamsplitter  to 
make  a reference  and  measurement  path.  The  energy  in  each 
path  is  measured  with  a separate  detector.  The  second  scheme 
uses  a reference  cell  and  a measurement  cell.  Each  cell  is 
alternately  placed  in  the  measurement  path.  Both  of  these 
schemes  are  currently  being  used  in  existing  instrumentation. 

The  beamsplitter  approach  is  being  used  in  a CO  monitor 
designed  for  NASA-Langley  by  TRW  to  remotely  measure  the  at- 
mospheric CO  concentrat ion  from  an  airplane.  This  system 
has  the  advantage  of  no  moving  parts  and  that  the  comparison 
between  the  two  beams  is  made  at  the  same  time.  The  disad- 
vantage is  that  a sophisticated  amplifier  gain  and  balance 
control  are  required.  The  minimum  modulation  which  they 
report  is  .01%. 

The  second  system  is  a rotating  cell  containing  two 
sections,  a reference  and  a balance,  which  are  alternately 
rotated  into  the  optical  path.  This  system  has  the  advan- 
tage over  the  current  modulation  scheme  in  that  no  energy 
is  lost  in  the  beamsplitter  and  the  difference  in  the  two 
optical  paths  is  minimal.  This  scheme  is  used  in  several 
instruments  built  by  Ford  Aeronautics  for  EPA  and  based 
on  the  data  collected  on  one  of  these  instruments.  The 
drift  in  M is  less  than  5%  of  the  expected  modulation  from 
N^O.  Hence,  we  are  reasonable  assured  that  this  scheme 
will  provide  satisfactory  results.  The  disadvantage  of  this 
scheme  will  provide  satisfactory  results.  The  disadvantage 
of  this  scheme  is  that  the  rotating  chopper  would  also  have 
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to  be  provided  for  the  HDO  and  CO^  measurements.  Based  on 
our  own  experience  and  on  information  received  from  NASA- 
Langley , there  is  a reasonable  expectation  that  with  an 
improved  chopper,  the  present  chopper  modulation  scheme 
will  give  satisfactory  results. 
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SUMMARY  AND  CONCLUSIONS 


| 

With  the  work  reported  here,  the  goal  of  developing 
a precise  understanding  of  atmospheric  molecular  absorp- 
tion of  DF  laser  radiation  has  been  attained  for  sea 
I level  paths.  Some  of  the  more  important  milestones 

* reached  during  this  two-year  effort  are  the  following: 

(1)  High  resolution  (s  .05  cm  *)  survey  spectra 
of  HDO,  N^O,  and  CH^  have  been  measured  and 

! compared  with  HI -TRAN  calculations  to  identify 

inaccuracies  or  omissions  in  the  AFGL  data  tape 
which  should  be  important  to  molecular  absorp- 
tion of  DF  laser  radiation.  These  spectra 
| identified  HDO  and  CH^  as  the  two  species  for 

which  the  molecular  line  data  base  was  inade- 
I quate  for  this  application. 

(2)  Synthetic  spectra  were  generated  at  the  long 

I wavelength  end  and  the  short  wavelength  end 

‘ of  the  low  power  laser  and  DDL  spectral  region  to 

, quantify  the  important  molecular  absorbers 

j 

| for  the  high  power  DF  devices  (NACL  and  MIRACL) 

and  for  broad  band  E-0  applications,  respectively. 

I (3)  A complete  base  of  high  resolution  HDO  measure- 

ments suitable  for  extracting  HDO  line  strength, 

I air  broadened  width,  and  position  parameters 

has  been  established. 

[ (4)  Algorithms  for  predicting  DF  laser  transmission, 

function  of  temperature  and  absolute  humidity 
| have  been  developed  for  sea  level  propagation 

paths . 

II 


i 
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(5)  HI-TRAN  modeling  for  predicting  infinite 
and  degraded  resolution  SMI  spectra  has  been 
developed . 

(6)  Modeling  of  SMI  spectra  and  DF  laser  trans- 
mission algorithms  have  been  verified  by  NRL 
field  data. 

(7)  HDO  content  in  the  atmosphere,  relative  to 
the  most  abundant  isotopic  variation  has  been 
found  to  be  typically  less  than  the  heretofore 
assumed  value  of  .03%  HDO/HHO.  This  has  been 
determined  by  NRL  using  a HDO  GFCS  designed 
and  constructed  by  SAI . 

(8)  A design  study  for  modifying  the  GFCS  to 
monitor  up  to  5 gases  important  to  DF  laser 
transmission  has  been  completed. 

The  work  reported  here  and  in  the  preceding  DF  work 
has  not  directly  addressed  aerosol  and  continuum  extinc- 
tion; nevertheless,  a better  understanding  of  these 
issues  have  been  achieved  since  the  modeling  developed 
has  allowed  continuum  extinction  to  be  extracted  from 
field  data.  Thus,  a powerful  tool  for  "backing  out" 
aerosol  and  molecular  continuum  data,  and  for  achieving 
predictive  models  for  this  extinction,  has  been  developed. 
At  the  present  time,  however,  it  must  be  concluded  that 
a verified  predictive  modeling  of  aerosol  extinction 
has  not  been  achieved. 

The  verification  of  the  variability  of  HDO  content 
as  a function  of  the  environmental  condition  (item  7)  has 
identified  the  need  for  determining  actual  atmospheric 
number  densities  of  the  ambient  as  well  as  the  variable 
gases,  rather  than  relying  on  the  accepted  global  average 
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values.  Construction  of  the  multi-channel  GFCS  (item  8) 
would  permit  monitoring  all  gases  that  are  most  relevant 
to  DF  HEL  test  and  operational  performance. 

In  addition  to  determining  the  actual  atmospheric 
molecular  content,  short  range  and  above  sea  level 
transmission  remains  as  the  unsolved  linear  propaga- 
tion issues.  Algorithms  for  molecular  absorption  for 
these  scenarios  have  not  been  developed,  nor  have  they 
been  developed  or  verified  even  for  sea  level  aerosol 
extinction. 
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APPENDIX 
HDO  DATA  LOG 


A log  of  all  181  data  sets  obtained  during  the 
present  program  has  been  maintained.  It  is  repro- 
duced here  as  an  Appendix  to  document  the  experimental 
conditions  and  other  information  relevant  to  its  use. 
The  log  is  described  in  Section  2. 
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